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FIELD OF THE INVENTION 

The present Invention relates to a method for estimating a propagation channel In the 
presence of transmit beamforming. It further relates to an estimator Implementing said 
5 method. 

Said method is suited to be used in particular within the UMTS standard and 
especially within a conventional RAKE receiver. 

BACKGROUND OF THE INVENTION 

10 A UMTS communication system, defined in the UMTS standard (see the hereby 

Incorporated reference "3GPP Technical Specifications, Rel. 99, 

http;//w WW .3qpp.orq/spe f -s^ p P rs, htmn , comprises at least a base station (BS) and a mobile 
phone terminal, referred to as user equipment (UE) throughout the document Both the BS 
and the UE includes a transmitter and a receiver. The BS sends signals by mean of the BS 
15 transmitter to the UE receiver in the downlink. The UE sends signals by mean of the UE 
transmitter to the BS receiver in the uplink. In uplink and downlink communication both the 
UE and BS map the data to be transmitted into logical channels. 

In particular the UMTS standard provides for the downlink communication two logical 
channels, a first channel denoted common pilot channel (CPICH) consisting of a continuous 
sequence of known pilot symbols and a second logical channel denoted dedicated physical 
channel '(DPCH). Said second logical channel DPCH consists of a first logical sub-channel 
denoted dedicated physical data channel (DPDCH), that carries the Information data symbols 
intended for the user of interest time multiplexed with a second logical sub-channel denoted 
dedicated physical control channel {DPCCH), that carries known user dedicated pilot symbols. 

Transmit beamforming can be employed at the BS during downlink communication to 
improve the reception at the UE of the user information data mapped Into the logical channel 
DPDCH. To perform transmit beamforming the BS needs to be equipped with multiple- 
antennas, in this way It is possible to provide optimal spatial filtering by shaping beams in the 
directions of the user of Interest while setting nulls In other users' directions to mitigate the 
interference to other users. Hence the use of transmit beamforming can significantly improve 
the efficiency In the spectrum usage, allowing a larger number of users to share the same 
spectrum, thanks to an improved user dedicated data reception and a reduced interference. 
While transmit beamforming is to be employed for the transmission of the logical channel 
DPCH it is not suited for the logical channel CPICH transmission that is to be uniformly 
broadcast in time and space (i.e. in all directions) to all usere. As seen from the UE receiver, 
the transmission of the logical channel DPCH occurs through a so called DPCH propagation 
channel consisting of the cascade of the BS transmitter chain, including the beamforming 
antenna system, the actual propagation channel from the BS to the UE and the UE receiver 
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chain. Unlike the logical channel CPICH downlink transmission occurs through a so called 
CPICH propagation channel consisting of the cascade of the BS transmitter chain. Including 
an omnidirectional antenna system In place of the beamformlng antenna system, of the 
actual propagation channel from the BS to the UE and the receiver chain. For the sake of 

5 brevity we shall refer to as DPCH channel and CPICH channel the above- DPCH propagation 
channel and CPICH propagation channel respectively throughout the present document. The 
UE receiver needs to estimate the DPCH channel in order to demodulate, decode, and recover 
the information data contained in the DPDCH. In the absence of transmit beamforming the 
CPICH channel and the DPCH channel are the same (up to a unlnfluentlal real scalar factor 

10 due to the different transmit power allocated to the DPDCH, DPCCH and CPICH logical 

channels). The continuous provision of known pilot symbols of the CPICH logical channel can 
be easily exploited at the UE receiver to estimate the CPICH channel, that is the same as the 
DPCH channel, and then use that estimate for DPDCH information data recovery. 
Unfortunately when transmit beamforming Is applied, the known pilots of the CPICH logical 

15 channel yet provide estimates for the CPICH channel but those estimates cannot be directly 
applied as channel estimates for the DPCH channel that being different from the CPICH 
channel. Then the UE receiver has to rely on different techniques to estimate the DPDCH 
channel in order to recover the information data. 

20 SUMMARY OF THE INVENTION 

Accordingly, it is an object of the invention to provide a method for DPDCH channel 
estimation in the presence of transmit beamforming optimally exploiting all information 
available from the logical channels DPCH and CPICH, while limiting computational complexity, 
and an associated estimator implementing said method. 

25 

To this end, according to a first object of the Invention, there is provided a method for 
estimating a propagation channel in the presence of transmit beamforming as claimed In 
claim 1. 

In addition, according to a second object of the Invention, It Is provided an estimator 
30 for estimating a propagation channel in the presence of transmit beamforming as claimed in 
claim 15. 

In a first non-limitative embodiment of the invention, the method is characterized as 
claimed In claim 2. 

In a second non-limitative embodiment of the Invention, the method Is characterized 
35 as claimed In claim 3. 

In a first non-limitative embodiment of the invention, the estimator is characterized 
as claimed in claim 16. 




In a second non-limitative embodiment of the invention, the estimator Is 
characterized as claimed in claim 17. 

As we will see In detail further on, this method enables the UE receiver to cope with 
5 transmit beamforming applied to DPCH transmission by jointly exploiting the information 
available at the UE from the CPICH and the DPCH logical channels. The first embodiment of 
the method optimally exploits that Information according to the maximum-a-posteriori 
criterion, which is well known to the person skilled in the art, to estimate the DPDCH 
propagation channel multipath coefficients as described above. We shall show that in the 

10 presence of transmit beamforming, each multipath component coefficient of the DPCH 

channel differs from the corresponding multipath component coefficient of the CPICH channel 
by an unknown complex scale factor already referred to as beamforming complex factor. The 
second embodiment of the method provides an estimate of that beamforming complex factor. 
Then an estimate of the multipath component coefficient of the DPDCH propagation channel 

15 is built by multiplying each multipath component coefficient estimate of the CPICH channel by 
the related complex scale factor estimate. In this way a considerable computation saving is 
possible with respect to the first method. Moreover the use of the CPICH channel estimates 
also allows to continuously track the channels variations with time due to the Doppler spread 
generated by the motion of the UE, thanks to the continuous provision of known pilot 

20 symbols of said CPICH logical channel. 



BRIEF DESCRIPTION OF THE DRAWINGS 

The invention and additional features, which may be optionally used to Implement the 
invention to advantage, are apparent from and elucidated with reference to the drawings 
25 described hereinafter. 



- Rg. 1 represents a diagram illustrating the main steps of said first method for DPDCH 
channel estimation according to the invention, 

- Fig. 2 represents a diagram illustrating the main steps of said second method for DPDCH 
30 channel estimation according to the Invention, 

- Rg. 3 illustrates an architecture schematic implementing the first method for DPDCH 
channel estimation according to the invention of Rg. 1, and 

- Fig. 4 illustrates an architecture schematic Implementing the second method for DPDCH 
channel estimation according to the Invention of Fig. 2. 
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DETAILED DESCRIPTION OF THE EMBODIMENTS 

1 Introduction 

In the following description, well-known functions or constructions by a person of 
ordinary skilled In the art are not described In detail since they would obscure the Invention in 
unnecessary detail. 

The present invention relates to a method and embodiments for the estimation of the 
propagation channel in a communication system when information data are transmitted 
through a beamforming system. Said methods are particularly suitable for Implementation In 
a RAKE receiver within the UMTS standard when information data are transmitted from a 
base station BS to a user equipment UE, here a mobile phone. 

The UMTS standard provides for the downlink communication two logical channels: a first 
channel denoted common pilot channel (CPICH) consisting of a continuous sequence of 
known pilot symbols and a second logical channel denoted dedicated physical channel 
(DPCH). Said second logical channel DPCH consists of a first logical sub-channel denoted 
dedicated physical data channel (DPDCH), that carries the information data intended for the 
user of interest, time- multiplexed with a second logical sub-channel denoted dedicated 
physical control channel (DPCCH), that carries known user dedicated pilot symbols. Every 
time slot a DPCH Is allocated, namely a DPDCH data field followed by DPCCH pilots field Is 
allocated. 

Transmit beamforming can be employed at the base station during downlink 
communication to improve the reception at the UE if the user information data mapped in the 
logical channel DPDCH. To perform transmit beamforming the BS needs to be equipped with 
multiple-antennas. In this way It is possible to provide optimal spatial filtering by shaping 
beams In the directions of the user of Interest while setting nulls in other users' directions to 
mitigate the interference to other users. Hence the use of transmit beamforming can 
significantly improve the efficiency In the spectrum usage, allowing a larger number of users 
to share the same spectrum, thanks to an improved user dedicated data reception and a 
reduced interference. While transmit beamforming Is to be employed for the transmission of 
the logical channel DPCH It is not suited for the logical channel CPICH transmission that is to 
be uniformly broadcast In time and space (i.e. In all directions) to all users. As seen from the 
UE receiver, the transmission of the logical channel DPCH occurs through a so called DPCH 
channel consisting of the cascade of the BS transmitter chain, including the beamforming 
antenna system, the actual propagation channel from the BS to the UE and the UE receiver 
chain. Unlike, the logical channel CPICH is transmitted through a so called CPICH channel 
consisting of the cascade of the BS transmitter chain, including an omnidirectional antenna 
system in place of the beamforming antenna system, of the actual propagation channel from 




the BS to the UE and the receiver chain. The UE receiver needs to estimate the DPCH channel 
in order to demodulate, decode, and recover the information data contained in the DPDCH. 
In the absence of transmit beamforming the CPICH channel and the DPCH channel are the 
same (up to a real scalar factor due to the different transmit power allocated to the DPDCH, 
5 DPCCH and CPICH logical channels). The continuous provision of known pilot symbols of the 
CPICH logical channel can be easily exploited at the UE receiver to estimate the CPICH 
channel, that is the same as the DPCH channel, and then use that estimate for DPDCH 
information data recovery. 

Unfortunately when transmit beamforming Is applied, the known pilots of the CPICH logical 
10 channel yet provide estimates for the CPICH channel but those estimates cannot be directly 
applied as channel estimates for the DPCH channel that being different from the CPICH 
channel. 

2 Fundamentals on Data Transmission and Reception in DS-CDMA 

In order for the information dedicated to a certain user to be transmitted the source 

15 Information data are properly mapped into a binary alphabet, i.e. Into bits, via source- 
encoding techniques. Then those bits are channel-encoded by using suitable channeJ codes in 
order to protect the information bits from all Impairment causes (noise, interference) during 
the transmission. The channel-coded bits are then modulated by mapping bits into symbols 
according to a certain modulation scheme (e.g. QPSK modulation, well know by the person 

20 skilled in the art) where each symbol Is represented by one or more bits (in QPSK modulation 
each symbol is represented by two bits). In DS-CDMA systems said symbols are then spread 
over a larger bandwidth by filtering with specific spreading sequences. In the literature 
spread symbols are referred to as chips Notice that the duration of a chip period is smaller 
than the duration of a symbol period by a factor approximately equal to the bandwidth 

25 expansion factor of the spreading operation. Said factor is equal to the number of chips per 
symbol period and it is referred to as spreading factor. Symbol sequences carrying 
information bits from different sources are spread with different spreading sequences in order 
for several users to share the same bandwidth and for the receivers to be able to discriminate 
and recover the desired user's symbol sequences based on the knowledge of the user specific 

30 spreading sequence. The obtained chips are then filtered by a pulse-shape filter (a root- 
raised cosine filter in the case of UMTS) and digital-to-analog (D/A) converted. The resulting 
analog signal Is then modulated to radio frequency and transmitted by the BS antenna(s). 
Similarly when transmit beamforming Is employed the same sequence of chips filtered by the 
pulse-shape filter is sent to all the antennas after multiplication by an antenna specific 

35 complex weight factor. The antenna specific weight factors can be collected In a so-called 
beamforming weight vector. The choice of said beamforming weight vector determines the 
shape of the Irradiation diagram, In terms of magnitude and phase, of the multiple-antenna 
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system of the BS. The signals transmitted from the antennas undergo several scattering, 
diffraction, and reflection phenomena through the actual propagation channel from the BS 
antennas to the UE receiver antenna. Noise and Interference from other users and 
neighboring BS's further contribute to Impair the useful slgnal(s) for the UE of Interest The 
5 radio frequency signal received at the UE antenna is demodulated down to baseband (or at 
Intermediate frequency IF) and analog-to-digital (A/D) converted to generate a digital 
baseband signal. Then the UE receiver processes that baseband signal In order to recover the 
useful information data intended for the user of interest For this purpose the receiver needs 
to estimate the cascade of the transmitter chain, the actual wireless propagation channel and 
10 receiver chain up to the A/D conversion, associated with the transmission of the useful data. 
In the literature that cascade is in general simply referred to as the channel. In the sequel a 
detailed rpathematical description of the above operations is provided. 

3 Signal and Channel Model 

General DS-CDMA signal and channel models are provided. Such models will then applied to 
15 the specific cases of UMTS CPICH and DPCH signals and channels. 

3.1 General Signal Model 

We consider a general case of transmission of a signal s(t) through a multipath channel wit h 
impulse response . The continuous-time complex base-band received signal prior A/D 
conversion at the receiver is modeled as 
20 rtO-xCO+vCO (1) 

where x(fy represents the part of received signal that comprises the useful data and v(t) 
denotes the noise plus interference term- The signal x{t) Is given by 



*(0= J/KM-r)z(r)dr ( 2 ) 

25 

where h(t,r) represents a time-varying channel impulse response and 

represents the transmitted signal after D/A conversion at the transmitter, withr c denoting 
the chip period and a(k) = s([k/M])d(k) , where M denotes the spreading factor, andjj 
30 denotes the flooring operator. Moreover s[n) represents the n-th modulation symbol (e.g., 
QPSK), d(k) denotes the /r-th chip of the spreading sequence, and yKf) represents the puise- 
shaplng filter that limits the system bandwidth. Note that this model can be used by systems 
with several spreading layers, like DS-CDMA system defined by the standard UMTS, IMT-2000 
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and IS-95 (see the hereby Incorporated reference "IMT-2000 htto://wwuy i m t-- 
ZOOO.orn/pnrh al/lndPxas p^ 

The channel h(f,r) Is assumed to follow the wide-sense stationary uncorrected 
scattering model with Rayleigh fading and multlpath response as described . hereby In the 
Incorporated referenced. G. Proakls Digital Communications, NY: McGraw-Hill 3rd edition 
1995" ' ' 

where /»denotes the number of multlpath components, and c/0 and r„ denotes the time- 
varying complex coefficient and the propagation delay associated with^-th multlpath 
component. The path complex coefficients c^t) vary depending on the velocity of the UE 
and of the surrounding scattering objects (e.g., buildings, hills but also other moving objects) 
with respect to the BS. When the BS uses Q antennas for transmit beamforming, the channel 
model expressed by (4) generalizes to 

Kt, r) = Z w^a^ ) Cp (t)5(j-T p ) (5) 

15 where 

" w = K w e _, f denotes the transmit beamforming weight vector 

a«?)=[«oW...,a Q .,(<9)] T denotes the antennas response vector in a direction*. 
The superscripts ( . )^and ( . )" denotes transpose and HermlHan transpose respectively. 
Since in a practical implementation the vector w typically may vary at most at the slot rate 
for the sake of notation Its variations with time are not explicitly taken into account within 
equation (5). 

By Inspecting equations (4) and (5) one may notice that, during a time slot, transmit 
beamforming only affects the amplitude of the path time-varying complex coefficients by 
a complex factor w"a(*„) and not the delays v in other words the channel's temporal 
structure Is not affected by the antenna response vector, neither by the transmit 
beamforming. Thus the path time-varying complex coefficient/^/)- W "a( W /)and CpW 
are the same in both case (4) and (5) up to a complex factor/?, = w »a ( * P ) representing the 
previously mentioned beamforming complex factor. Therefore a channel estimate according 
to (5) can be built from estimates of the delays r„of the coefficients and of Afcr jp-i, .... 
30 P, 
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3.2 Received Signal Model for DPCH and CPICH 

We now apply the previous general channel model to the DPCH and CPICH channels 
considering the case of UMTS where the CPICH and DPCH channels are superposed. Then 
the continuous-time complex base-band received signal (1) can be written as 

where ,^(0 and are the signal components associated with the CPICH and DPCH 
channels respectively, given by 

JWM- 



where 



p-i 



and 
where 

20 - ^note the n-th modulation symbol associated with the CPICH and the Hh 

modulation symbol associated with the DPCH respectively, and d^k), d JjKh {k) donote 
the fc-th chips of the CPICH and DPCH spreading sequences respectively. 

3.3 Remarks 

It will be noted that the spreading factor M may be not the same for the CPICH and DPCH 
25 logical channels. In the UMTS standard the CPICH has a fixed spreading factor equal to 256 
whereas the DPCH has a spreading factor that can take on values in the set {4, 16, 32, 64, 
128, 256, 512}. 

Furthermore, the CPICH symbols are all the same, I.e., ^(n)=AH) for any «, 
where A is a real amplitude factor taking into account the power that the BS assigns to the 
30 CPICH logical channel. The a-priorl knowledge at the UE receiver of the symbols comprising 
the CPICH can be exploited to estimate the channel ft^ as expressed by (4). In the absence 
of transmit beamformlng on the DPCH k+> = ft** so that such an estimate based on the 
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CPICH symbols would suffice for the receiver to recover the desired user's Information 
symbols on the DPCH. On the contrary, If transmit beamformlng Is applied to the DPCH the a- 
priori knowledge of the CPICH symbols yet allows to easily estimate the channel h vleA as 
expressed by (4), but not the channel associated with the DPCH expressed by (5) which 
differs from the due to the presence of the unknown coefficients fip. Hence other 
approaches need to be addressed to provide an estimate of the channel h pch necessary to 
recover the information symbols transmitted over the DPDCH. 

We have implicitly assumed above a digital baseband receiver where baseband signals 
are processed after A/D conversion. In other words the signal expressed by (6) is A/D 
converted, prior all baseband processing, in particular we consider the UE receiver to be a 
discrete-time RAKE receiver (see the hereby incorporated reference W J. G. Proakls Digital 
Communications, NY: McGraw-Hill, 3rd edition, 1995"). We shall notice that the methods 
hereby addressed can be applied to all types of baseband receivers relying on channel 
estimation though. 

4 Channel Estimation 

Since the delaysr, are not affected by the presence of transmit beamforming we shall 
assume that they are already recovered by usual delay estimation/tracking techniques (e.g. 
delay lock loop, see the incorporated reference W J. G. Proakls Digital Communications, NY: 
McGraw-Hill, 3rd edition, 1995" exploiting the knowledge of the CPICH symbols) and we shall 
focus on the estimation of the channel path coefficients c„(/) and fi„. In order for those 
coefficients to be estimated we consider the signals at each finger of the RAKE receiver after 
despreading of both the CPICH and the DPCH. 

For the sake of simplicity we shall limit the treatment to the single path/finger case, 
the muiti-finger case being a straightforward generalization of It Thus, we will drop the 
subscripts to denote the time-varying complex coefficient c p (t) hereinafter also referred to as 
path fading coefficient, or simply path coefficient. Moreover for the sake of brevity in the 
sequel we shall refer to as channel estimation and channel estimate to as the channel path 
coefficient estimation and estimate respectively. 

Furthermore, in a first Instance, we also limit the analysis to the case where the 
DPCH spreading factor M Is equal to that of the CPICH. We shall generalize the analysis to 
other spreading factors successively. 

In addition we need to model the previously mentioned interference and noise term. 
In particular the multi-access Interference that arises from the Interference among different 
users and the interferences due to the CPICH over the DPCH and/or due to different data 
flows dedicated to the same user, and the Inter-path Interference that Is due to the 
Interference between the Impulse responses associated with paths at different delays r . 
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Assuming that the whole Interference and noise to be modeled as a zero-mean 
complex circularly symmetric additive white Gaussian noise (AWGN), the CPICH and DPCH 
de-spread signals, at a RAKE finger are respectively given by 

y^*(i»)-c(»)***»C»)+»w(") (7) 

where v^n) and v«m0O denote the noise-plus interference terms associated with the DPCH 
and CPICH respectively, assumed to be independent, c(n) is the discrete-time representation 
at symbol rate of the path time-varying complex coefficient previously denoted <*i) (or c p (0), 
,4*00 and Sepl c(n) denote the symbols of the DPCH and CPICH logical channels respectively, 
at a symbol time n. Under the assumption that the CPICH and DPCH channels may vary every 
symbol time «, and that the interference and the noise can be modeled as a zero mean 
AWGN, we then build a maximum likelihood (ML) static Instantaneous channel estimates at 
every symbol time » as follows 

^W-iUW^W-^+W") (8) 
= 'iaOOXo-W = w"a(0)c(n) + v^(») 

where, /(«) denotes the complex conjugated of *(»), v^(n) and V<M (») represent the 
AWGN associated with the CPICH and DPCH channels respectively, and where, without loss of 
generality, we assumed the symbols of both channels to be normalized such that |s^*(n)| = 
IW«)I = 1 • In order for the est" 113165 W to te computed one needs to know the 
transmitted symbols s' ipth (n) and <«(«) . The symbols are all known at the receiver. 
Within the DPCH only the dedicated pilot symbols j^OO comprising the DPCCH are known 
but not the data symbols comprising the DPDCH. Therefore if DPDCH symbols s^fr). are to 
be used in (8) they first need to be estimated by, e.g., a decision-direct technique. The ML 
estimates (8) are therefore formed from the observations y,** and y^ at the output of a 
finger in a data aided (DA) fashion when there is provision of a priori known pilot symbols, as 
for the cases of CPICH and DPCCH, and in a decision direct{DD) fashion In the presence of 
unknown data symbols as for the DPDCH. The DD approach requires the use of hard 
decisions taken on the DPDCH at the RAKE receiver output as pilot symbols. Therefore, from 
(8) we have 

WO - *l«kO0y+*<r0 = TY H *W)c(n) + v^(«) (9) 
for the DPCCH case (DA) and 

W»> = $U(»)jW») = w " a (*M»> + W»> (10) 
for the DPDCH case (DD) where S^iri) represents the estimate of the «-th DPDCH symbol 
built as a hard decision taken at the RAKE receiver output Equations (8)-(10) represent the 
basis for the derivation of the DPDCH channel estimation algorithms addressed in the sequel. 



In order to support the next derivation further details are provided concerning the physical 
behavior and related models for the CPICH and DPCH channels. 

4.1 Physical Interpretation of CPICH and DPCH Fading Channel Models 

Due to the large bandwidth of the UMTS W-CDMA system the channel from the BS to a UE 

5 can be modeled as comprising a set of discrete multipart! components. In particular it is 
possible to approximate the channel impulse response as a linear superposition of a finite 
number of independent multipath components as expressed by equation (4) and (5)- In the 
preceding we have also referred to a multipath component simply as a channel path. A path 
Is characterized by a delay and a complex coefficient When the UE (and/or any other 

10 scattering object) Is in motion with respect to the BS each path complex coefficient may vary 
with time. Thus the overall channel impulse response varies with time as well, as expressed 
by equation (4) and following. This phenomenon is typically referred to as Doppler spread, 
from the Doppler effect arising when in a communication either the transmitter or the 
receiver moves one with respect to another. The channel Doppler spread represents the rate 

15 at which the channel impulse response varies, or focusing on a channel path, the rate at 
which the path complex coefficient varies with time. The Doppler spread Increases with the 
velocity of the UE with respect to the BS. From a physical standpoint to each path of the 
CPICH channel corresponds a path on the DPCH channel with the same delay r and the 
same Doppler spread but with different complex paths' amplitudes, namely c(n) and fic(n) 

20 respectively, according to equations (4) and (5). Then as long as the factor^ is constant, both 
CPICH and DPCH channel paths vary at the rate of c(n) f namely they share the same Doppler 
spread. Hence the two complex paths' amplitudes temporal variations are the same while 
they always differ only by the beamformlng complex factor p that is determined by the scalar 
product of the beamformlng weight vector w and the antenna array response a(0) In 

25 correspondence of the direction 0 of that path. 

4.2 Optimal Joint MAP CPICH-DPCH Based DPDCH Channel Estimation 

The optimal joint CPICH-DPCH MAP DPDCH channel estimation problem Is formulated as a 
direct extension of the simpler CPICH MAP channel estimation problem addressed in the 
hereby incorporated reference W J. Bsltersee, G. Fock, P. Schultz-Rittich, and H. Meyr, 
30 "Performance analysis of phasor estimation algorithms for FDD-UMTS RAKE receiver," IEEE 
Symp. on Spread Spectrum Techonologtes and Applications, NJIT, NJ, September 2000." 
The general formulation of the MAP optimization criterion corresponds to the following 
optimization problem 

Cmap (*) = argmax p(c(k) | y , s) 
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namely c MP {k)\s the argument that maximizes the a posteriori probability density function 
P(c(*)|y,s)°f c(ft) given the vector yof the observations ^(«)and^(n) expressed by 
(7) and the vector s of the pilot symbols*^ (w) and (») . 

We shall see that even though optimal MAP DPDCH channel estimation is generally too 
computationally demanding for a feasible implementation within a realistic UE, its derivation 
provides the proper guidelines for other reduced complexity estimator for the DPDCH channel 
yet according to certain optimatity criteria. 

Under the assumption that the channel paths are Gaussian complex random process, the 
maximum a posteriori probability (MAP) DPDCH channel path coefficient estimate at time 
instant ft is given by the conditional expected value (ft) of e M (ft) = c apdeh (ft) given the 
observation column vector c a comprising all the available ML instantaneous channel path 
coefficients estimates c^(«) with m*ft,and B^Qi) defined according to equations (8)- 
(10), i.e., 

Equation (11) can be simply explained by observing that the ML estimates c M (ft) according to 
equations (8)-(10) are still Gaussian complex random variables with mean given by the true 
DPDCH channel (ft) = ftd/c) and variance given by the variance of the noise term (ft) . 
In order to achieve a MAP estimate the a-posterfori probability density function of c UL (K) (I.e., 
the probability density function of c M (k) conditioned to c ML ) has to be maximized. It is well 
known that the maximum of a Gaussian probability density function is reached in 
correspondence of the mean value of the distribution. Henceforth taking the expectation in 
equation (11) maximizes the conditional probability density function yielding a MAP channel 
estimate. It can be shown (see e.g. the incorporated referenced. Meyr, M. Moeneclay, and 
S. Fetchel, Digital Communication Receivers: Synchronization, channel estimation and signal 
processing, John Wiley & Sons, New York, 1998") that a MAP optimal channel estimate (11) 
can be built by low-pass filtering a set of Instantaneous ML channel estimates (8) as follows 

W*)= fW£ *«. (12) 
wheref=[/W-l) ... /(0)fln (12) Is the vector of the low pass filter coefficients, with 
N-K^+K**, *^and Kg** being the numbers of the available ML Instantaneous 
channel estimates of the DPCH, c^m), and of the CPICH, c^(»), respectively. Under the 
previous assumption that the channel paths are Raylelgh fading the optimal filter 
f = [f (N _ d /( 0 )f yielding a MAP channel estimate Is given by 



(13) 
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where R a = E^cj^} represents the covariance matrix of the ML channel estimates, 
V represents the vector of the cross-correlation between the ML channel 

estimates and the true DPDCH channel at time Instant k . 

The optimal filter f given by equation (13) corresponds to the well-known Wiener filter, 
5 also known as a linear prediction niter. Depending on whether the vector 6^ Includes only 
past or only future, or both past and future ML instantaneous channel estimates of the form 
(8)~(10) with respect to the DPDCH channel c^k) , equation (13) gives rise to a forward, 
a backward, or a forward-backward linear prediction filter respectively. 

The optimal linear prediction filter f can be easily computed also in the more general case 

10 of Rice fading channels as shown in deep detail in the hereby incorporated reference ™H. 
Meyr, M. Moeneclay, and S. Fetahel, Digital Communication Receivers: Synchronization, 
channel estimation and signal processing, John Wiley & Sons, New York, 1998." Despite the 
theory clearly indicates a closed form expression for the optimal solution to the problem of 
MAP channel estimation, as demonstrated above, a practical implementation of such an 

15 optimal estimator poses several complexity issues which generally prevent Its feasibility within 
a practical UE architecture. In the next section we shall address the Implementation of a MAP 
channel estimator. To this end we shall rely on the concepts of Wiener filtering, linear 
prediction, as well as adaptive and LMS-Kalman filtering oriented to MAP channel estimation 
along with fading channel models, which are considered to be well-known by the person of 

20 ordinary skill in the art. 

4.2.1 Implementation of the Optimal DPDCH MAP Channel Estimation 

In order to Implement the optimal Joint CPICH-DPCH MAP DPDCH channel estimator in a 

practical UE the optimum Wiener FIR filter coefficients f =r/(#-i) <m /(0)] r are to be 

25 computed. For that computation to be possible the statistics of the observation c ML and of the 
true channel c^^are needed at the UE. These quantities can be estimated assuming 
structured channel and signal models, such as the Jakes' fading model (see the hereby 
incorporated reference "J. G. Proakis Digital Communications, NY: McGraw-Hill, 3rd edition, 
1995"), and AWGN model for the noise and interference terms. Then estimates of the model 

30 parameters are needed, i.e., the Doppler spread (or the UE velocity relative to the BS), the 
signal power and the Interference-plus noise power spectral density. To save computations, 
the UE may select the optimal Wiener FIR filter from a set of pre-computed filter coefficients 
based upon the estimates of the parameters of the statistic channel model. Unfortunately 
even neglecting the complexity involved by the parameters estimation, this approach Is 

35 unfeasible. Indeed the number of filter coefficients to be stored would be overwhelming 
compared to the memory storage capabilities of a real UE. A sub-optimal solution would 
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consist of using a reduced set of Wiener filters that can cope with the maximum envisaged 
channel Doppler spread (normalized with respect to the symbol duration). Different filters 
would be selected based only upon the slgnal-to-noise-plus Interference ratio (SINR) estimate 
and no longer on the Doppler spread estimate. Indeed, estimating the Doppler spectrum is 
5 generally much harder than estimating the SINR. Moreover SINR estimates are needed 

anyhow to perform dosed loop power control. In order to guarantee that no Information loss 
occurs at the maximum Doppler spread for a given spreading factor, these filters need to be 
designed to cope with the widest envisaged Doppler bandwidth. This solution yields 
significantly suboptimal performances at low Doppler spread (i.e. low UE velocities) and/or 
10 low spreading factors, because the filter bandwidth results much larger than the one of the 
channel path process in these cases, leading to poor noise suppression. 

An alternative approach consists of building unstructured estimates of the observation 
covarlance matrix R*and the channel autocorrelation sequence v , without assuming a 
specific statistic model. Once such estimates are available the linear prediction filter f can be 
15 built according to equation (13) by replacing R a and r & . with their estimates. In practice ,t 
would be more convenient to estimate dlrectiy the Inverse of R, to avoid the complication of 
inverting the estimate of the R s . For this purposes several techniques are available in the 

A rTore practical solution when there Is provision of known pilots, I.e. over the DPCCH 
20 and/or CPICH, or reliable hard decisions taken via DD mechanism over the DPDCH are 

available, consists of estimating directly the Wiener filter coefficients in adaptive fashion, by 
sufficient training. Adaptive filtering inherently provides tracking of the short-term and long- 
term channel variations. Tracking and adaptation of the prediction filter can be performed m . 
a continuous fashion or on a biock basis, e.g. by adapting the FIR prediction filter every slot 
25 A lower complexity approach than FIR Wiener FIR linear filtering consists of using a first 
order HR LMS-Kalman as linear predictor. The IIR loop gain factor needs to be smaller than 
one for the IIR filter to be stable. Under the assumption of a parametric channel statistic 
model it can be shown that the optimal loop gain factor yet depends on the noise-plus 
interference term and the channel path statistics. Moreover, the optimal loop gain cannot be 
30 rendered smaller than one for high Doppler spread and large spreading factors, as it is 
' demonstrated In the Incorporated reference «3. Bsltersee, G. Fock, P. Schultz-Rlttich, and H 
Meyr "Performance analysis of phasor estimation algorithms for FDD-UMTS RAKE receiver, 
JEEEfSymp. on Spread Spectrum Teehonologles and Applications, NJIT, NJ, September 
2000 - Therefore, in order to be stable the LMS-Kalman linear predictor can no longer be 
35 optima, in those cases. Nevertheless the LMS-Kalman linear predictor is still an optimal 

solution for low UE mobility (Doppler spread smailer than 222 Hz) and small spreading factor 
(less than 128). 
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The above considerations allow estimating the amount of complexity involved by the 
computation of the optimal prediction filter. That complexity represents a major hurtle for a 
feasible implementation of the MAP channel estimator within a realistic UE so that other sub- 
optimal but less computationally intensive solutions need to be addressed. The presence of 
5 the expectation operator In the MAP estimator expressed by equation (11) indicates that any 
sub-optimal estimator should low-pass filter the instantaneous ML channel estimates 
expressed by equations (8)-(10). The design and implementation of a sub-optimal low-pass 
filter for the ML channel estimates (8)-(10) in place of the optimal Wiener filter has to 
achieve a trade-off between complexity saving and performance degradation with respect to 
10 the optimal MAP estimator. Very simple low-pass filter structure can be shown to achieve 

nearly optimal performances in most practical scenarios with significantly reduced complexity. 
These considerations are all about the practical design of a channel estimation mechanism 
suitable for a feasible UE. In the light of the previous derivation we shall Introduce a low 
complexity but still nearly optimal DPDCH channel estimation algorithm. 

15 4.3 Joint CPICH-DPCH Based Low Complexity DPDCH Channel Estimation 

In section 3.1 we have seen that the channel paths associated with the CPICH and DPCH 
share a common structure. In particular, it has been shown that both channels associated 
with the CPICH and DPCH behave in the same way up to said beamformlng complex factor p 
due to the transmit beamformlng. The knowledge of said complex factor £for all channel 

20 paths allows building an estimate of the DPDCH channel from an estimate of the CPICH 
channel. This way the DPDCH estimation problem splits In two simpler estimation problems, 
one of estimating £and another of estimating the CPICH channel c(k) . Optimal MAP channel 
estimation can be considered to estimate the CPICH channel c(k) with significant complexity 
reduction with respect to the complete joint CPICH-DPCH MAP DPDCH channel estimation 

25 problem. The optimal MAP CPICH channel estimate at time Instant k is given by the 
expectation of the ML Instantaneous CPICH estimate c^(ft) given all available ML 
Instantaneous CPICH estimates of the fbrmc^W for* * k . The solution to this problem has 
the same form as (12) and (13), where the observation column vector c M Is to be redefined 
so as to comprise just the CPICH and not the DPCH ML Instantaneous channel estimates 

30 (8)— (10), and the cross-correlation vector In equation (13) is to be redefined as 

V - BCc^c'C*)} . A detailed derivation of a MAP CPICH channel estimator Is provided in the 
hereby incorporated reference n J. Bsltersee, G. Fock, P. Schultz-Rittich, and H. Meyr, 
"Performance analysis of phasor estimation algorithms for FDD-UMTS RAKE receiver," IEEE 
6P Symp. on Spread Spectrum Techonologles and Applications, MIT, NJ, September 2000." 

35 Similar considerations as those mentioned In section 4.3 apply also In this case. If the UE 
cannot still afford the complexity of a MAP CPICH channel estimate lower complexity 
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approach can be pursued by adopting sub-optimal low pass filtering in place of the optimal 
linear prediction filter. In the sequel we shall refer to 2^(4) as an estimate of the CPICH 

channel c(k) regardless the low pass filter employed for the estimation. Along with an 
estimate (k) of the CPICH channel c(k) , the UE has to compute an estimate fi of said 
5 complex factor /?, to build an estimate of the DPDCH channel as 2^ (fc) « fic^ (k) . To this 
end, we shall derive an optimal ML estimator for/7 that can exploits both the DPCCH pilot 
symbols in DA fashion and/or the DPDCH data in a DD fashion, jointly with the pilot symbols 
provided by the CPICH according to equations (8)— (10). 

10 4.3.1 Maximum Likelihood Estimation of the Beamfbrmlng Complex Factor p 

In order to build a ML estimate for the factor/? we need to build ML estimates of the cross- 
correlation between DPCH and CPICH instantaneous ML channel estimates, and of the 
autocorrelation of the CPICH Instantaneous ML channel estimates. For Oils purpose we shall 
rely on the expressions of the ML instantaneous DPCH and CPICH channel estimates (8)— 
15 (10). The theoretical expressions for the DPCH-CPICH ML channel estimates cross-correlation 
and CPICH ML channel estimate autocorrelation are 

<U0 =B{c^(»)5; Wj (n~0} (M) 

and 

* m (f) - Eic^Qt^in -/)} (15) 

20 respectively, where a correlation lag / * Ois to be selected for noise suppression, as it will be 

shown In the sequel. From equation (8), we have 

c**dX*)«c(«)+ v q*ft00 

c dpcb (n) = w"a(0)c(«) + v dpch (*) 

where we denote p - w K a(0) . Under the assumption that the noise terms (») and 
v ^(n)are independent zero-mean additive white Gaussian noises, expressions (14) and (15) 
25 for/*oyleldto 

4(0-B^(»)c; w (n-/)} = 

« J3B{c(n)c (n - /)} + >SE {c(n)v^(« - 0} + E{v^(«)c\»-/)}+B{v rf/KA (n)v;^(n-/)} = 
==/fc{c(«)c>-/)} 

and 

= E{c(n)c (*-/)} 
respectively. Thus the complex factor £ Is readily given by 
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for any/*o (16) 
A practical Implementation requires the expectation operator E{.} to be replaced by a feasible 
estimator yielding estimates^/) and 4(/)for the cross- and auto-correlation (14) and (15) 
respectively, to build an estimate/? for /? according to equation (16). A simple and effective 
estimate for^(/) and ^(/)is given by the sample temporal average, i.e. 

0*(O=-j= — T. a *<rOK r kk(«-Q (17) 

and 

&c(0 = — Z a *«( B ) s i*»(»-0 (18) 

where n = m denotes the time Instant which the averages (17) and (18) start from. Notice 
that for / > 0 it Is needed s +l . The estimates (17) and (18) are consistent 
estimates and commonly known as sample estimators. They are also known to be ML 
estimates of the cross and auto-correlation (14) and (15) under the assumption of Gaussian 
observations, as In our case (see the hereby Incorporate reference W B. Porat, Digital 
processing of random signals, Englewood Cliffs, NJ, Prentice-Hall, 1994"). They are also 
known to be biased estimates for finite ,-and asymptotically unbiased, i.e. unbiased 
as*** goes to infinity. Moreover It is well known that the ratio of two ML estimates is still a 
ML estimate. Thus a ML estimate for J3 Is given by 

- , A -j?— 2£**(")c^w0i-Q 

A 1 ) "-dpch n-m 

At) 1 (19J 



K 



cplch 



where the dependence on the lag / is due to the finite number of symbols tf^and 
^available In practice. Notice that In the case where = jc^ the estimator (19) takes 
the form 



X m Z^(")c^*( M -0 



^ (/) = ^. = _^_ ( 20) 



The estimator (20) has the advantage of avoiding the need of computing the ratio 
K cpich I K dpch and therefore has a reduced complexity with respect to the general estimator 
(19). In order to further limit the complexity of estimator (20) one might avoid the presence 
of a DD mechanism to build the ML Instantaneous channel estimates over the DPDCH 
according to equation (10). In this case only the known pilot symbols comprising the DPCCH 
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would be used to form DPCCH ML Instantaneous channel estimates expressed by equation (9) 
and it would result ■» K dpixH . 

Since equation (16) holds for any /*0, A.© expressed by (19) and (20) represents 
a ML estimate of p for any I * 0 . It Is also well known that a linear combination of ML 
estimates is still a ML estimate. Thus in a preferred non-limitative embodiment, in order to 
improve the estimate of fl , one may consider the weighted sum of ML estimates of the 
cross- and auto-correlation at different correlation lag I , i.e. 



A _ l=\ M n~m {Z±) 

PhiL ~ L L 



where the coefficients a, and b, are to be chosen according to the accuracy 
10 (expressed in terms of mean square error) of each ML estimates ^(/)and MO (see the 
incorporated referenced. Porat, Digital processing of random signals, Englewood Cliffs, NJ, 

Prentice-Hall, 1994"). 

Moreover it can be shown that the estimators (17) and (18) are asymptotically 
efflclentonM for the first correlation lags /=1,...,I, where I depends onthe statistical 
15 properties of the random processes c^(»)and c^(») (namely for l>L the variance in the 
estimate even asymptotically does not attain the Cramer-Rao lower bound, a well known 
concept to the person skilled in the art, as it is shown hereby in the Incorporated reference 
«B. Porat, Digital processing of random signals, Englewood Cliffs, NJ, Prentice-Hall, 1994"). 
In another preferred non-limitative embodiment, the following weighted sum of ML 
20 estimates can be considered 

A.-±riA.« (22) 

M 

where a: represents the maximum correlation lag at which the cross- and autocorrelation 
(17) and (18) are estimated. Similarly to the estimator (21) the optimal weighting coefficients 
Yl can be computed based on the accuracy (expressed in terms of mean square error) of 
25 each ML estimates 0 m 0). in practice y, will increase with the accuracy of the estimate 
p M (f), rendering the contribution of the most accurate estimates to the sum (21) more 
significant than that of less accurate estimates. Similarly to estimator (21) the efficiency of 
the estimates^/) will decrease as the lag I Increases. Thus also in this estimator a 
maximum lag £ can be determined for the estimate &*(/) to be asymptotically efficient 
30 (see the hereby Incorporated reference »B. Porat, Digital processing of random signals, 
Englewood Cliffs, to, Prentice-Hall, 1994") where we remark that it would be In general 
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K * L . The derivation of the optimal weight coefficients a, and b, , or y, although of 
theoretical Interest, results far too complex to allow a simple dosed form expression suitable 
for a realistic implementation hence beyond the scope of this Invention, However it can be 
shown that since the accuracy of the estimates p m (f) as well as the accuracy of the 
5 estimates $ de {i) and generally decreases as the correlation lag / increases the 
coefficients y l as well as a, and b t will decrease accordingly. 

In the light of the above considerations the least complex ML estimator with the 
maximum estimate accuracy Is that of equation (20) with / = 1 , i.e. 

'*®"Im"v* : (23) 

n«ro 

10 Notice that estimator (23) has also the advantage of minimizing the delay needed to compute 
the cross and the autocorrelation at numerator and denominator respectively, reducing the 
complexity and the need of additional memory storage. 



4.3.2 Practical Means for Improving and Tracking the Beam forming Factor Estimate 

15 There might be practical circumstances where the beamforming factor can stay constant 

over several slot periods. That may occur when the BS optimizes the beamforming weight 
vector w on the basis of the paths angles O p and the paths' average energies. In such cases 
since those quantities can stay approximately constant over several slot periods the BS will 
not modify the beamforming weight vector w so that the factor p p = w^a^) for 

20 p will stay approximately constant during the same time period. That fact can be 
exploited to improve the accuracy of the estimates (19)— (23) by averaging/filtering the 
estimates of the cross- and auto-correlation at fixed lag / * 0 (17) and (18) computed over 
successive slot periods/ namely over the period where the factor p is constant. Then for that 
averaging-filtering operation to be properly implemented the UE has to be informed 

25 whenever the beamforming factor £ results abruptly modified either because the paths' 

angles have changed and/or the BS has modified the beamforming weight vector w , in order 
for the UE to reset/ reinitialize the averaging-filters memory. Indeed if the averaging/filter 
memory is not reset the UE would account for cross-correlation and autocorrelation estimates 
(16)— (17) corresponding to an out-of-date beamforming factor/? . However one may notice 

30 that there might be no need of resetting the averaging-filters memory in the case of slow 
variations of the beamforming factor p . Indeed appropriate filtering of successive cross- and 

auto-correlation estimates (17) and (18) inherently performs the tracking of slow variations of 
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the beamforming factor p . For this purpose we suggest the implementation of a low 
complexity first-order low pass IIR filter with a suited forgetting factor (first-order IIR 
filtering is well known by the person ordinary skilled in the art). Eventually in order for the 
BS to notify the UE about abrupt changes in the transmit beamforming a signaling overhead 

5 is necessary though- If the standard under consideration, e.g. UMTS, does not envisage to 
support such signaling from the BS to the UE, either the averaging filters have to be removed 
or the UE needs to implement an additional mechanism able to detect whether an abrupt 
change occurred in order to make a decision on whether the averaging filter memory is to be 
reset Such an additional detection mechanism as well as the cross- and auto-correlation 

10 averaging-filters, clearly come at the expense of complexity increase that might Impact the 
UE design. However if due to complexity limitation the averaging-filters are not Implemented 
the BS will not have to signal abrupt beamforming factor changes to the UE avoiding 
additional network overhead. The UE shall then up-date the estimate of the beamforming 
factor every slot, independently by the previous slots estimates. Furthermore one may 

15 observe that in this case the UE can keep estimating the factor fi even independently of the 
presence of transmit beamforming. Indeed when no transmit beamforming is performed at 
tfte BS, fi will reduce to a real factor representing the transmit amplitude offset between the 
CPICH and DPCH. 



20 43.3 DPDCH Channel Estimation at Different DPCH Symbol Rates 

So far, for the sake of simplicity in the derivation we have assumed as a particular instance 
the spreading factor of the DPCH to be equal to 256, in order for the DPCH and the CPIGH to 
have the same symbol rate (i.e. the same number of chips per symbol). However as 
previously mentioned the spreading factor of the DPCH can take on values in the set {2, 4, 8, 
25 16, 32, 64, 128, 256, 512} so as the DPCH symbol rate generally differs from that of the 
CPICH whose spreading factor is always fixed equal to 256. Hence in a preferred non 
limitative-embodiment, the instantaneous ML channel estimates (8)— (10) generalize as 
follows 

= » c(n) + v epleb (n) 

30 where, given the DPCH spreading factor M , the index / corresponds to the / -th DPCH 

symbol spread over M chips, and n corresponds to the n -th CPICH symbol always spread 
over 256 chips. Then since to any CPICH symbol correspond 256/ M DPCH symbols In order 
for the previous algorithm to be used either the CPICH or the DPCH ML channel have to be 
interpolated to match the symbol rate of the other channel estimates. 

35 Moreover In deriving the ML instantaneous channel estimates (8)— (10), both the 

CPICH and the DPCH channels are implicitly assumed constant over the longer symbol period 
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of the two said channels, namely the Doppler spread Is always assumed significantly lower 
than the lowest symbol rate allowed by the standard of interest We shall briefly show that 
this assumption holds in practice by considering. the case of UMTS standard. Indeed the 
lowest symbol rate arises when the spreading factor Is maximal, i.e. for M =512 . In the 
UMTS standard the chip rate is fixed equal to 3.84MHz (see the Incorporated reference tt 36PP 
Technical Specifications, Rel. 99, http://www.3oDP.ora/<;ppr«;/spe f ^.htmi^ r so that the lowest 
symbol rate for the DPCH Is equal to 3.84MHz/512= 7.5 kHz. The UMTS standard also 
envisages supporting wireless mobile services where the UE velocity can be up to 500km/h, 
l.e. In the presence of a Doppler spread at most equal to 1kHz that Is still largely below 
7.5kHz for the above assumption to be valid. Then in a preferred non-limitative embodiment, 
interpolation Is done to adapt the ML Instantaneous CPICH and DPCH estimates to the lowest 
symbol rate between said CPICH and DPCH. In particular, when the DPCH spreading: 
factor M Is smaller than 256, the interpolation to the CPICH symbol rate is obtained by 
averaging DPCH ML static channel estimates over a number of symbols equal to 256/M, 
whereas when M=512, being the DPCH symbol rate lower than that of the CPICH, the 
average Is taken on the CPICH over two consecutive ML Instantaneous channel estimates. 
One can notice that this embodiment have the advantage of not coloring the estimation 
noise, I.e. noise samples after Interpolation are still uncorrected so that equation (16) still 
holds. 

4.4 DPDCH Channel Estimation Algorithms 

In this section we summarize the steps of the joint CPICH-DPCH MAP DPDCH channel 
estimation algorithm as well of the joint CPICH-DPCH low complexity DPDCH channel 
estimation algorithm derived In sections 4.2 and 4.3 respectively. 

4.4.1 Optimal Joint DPCH-CPICH MAP DPDCH Channel Estimation Algorithm 
In this section we provide a summary of the main steps to Implement the optima joint DPCH- 
CPICH MAP DPDCH channel estimation algorithm described above. The following steps are to 
be performed to provide an optimal MAP DPDCH estimate according to (11). 

1. Compute the ML instantaneous CPICH and DPCH channel estimates («) and 
2^(0 according to equations (8)— (10) and the generalization (23). 

2. Perform interpolation of the above obtained ML Instantaneous DPCH and CPICH 
channel estimates to the lowest symbol rate of said DPCH and CPICH, i.e. 

• if the DPCH spreading factor M < 256 the DPCH ML channel estimates are to be 
Interpolated in order to match the CPICH symbol rate 
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. if the DPCH spreading factor M - 512 the CPICH ML channel estimates are to be 

Interpolated in order to match the DPCH symbol rate 
• if the DPCH spreading factor M = 256 no interpolation is needed 
The so interpolated ML instantaneous channel estimates are collected in the following 
vectors 

To further reduce the complexity we limit the estimates c^n-K^ +l),....c^(»)to 
the estimates over the DPCCH c^i-K^ +l),..,a* crt (0 in order to avoid the 
implementation of a DD mechanism necessary to provide pilot symbol estimates over 
the DPDCH. The vector =[c^ c^f Is formed. 

3. Build the optimal linear prediction filterf according to equation (13) either by 
implementing one of techniques described in section 4.3 

4. Filter the ML instantaneous channel estimates c m =[c^ c^] r obtained at step 3 
with the linear prediction filter computed at step 4 according to equation (12) to 
obtain a MAP DPDCH channel estimate c apdch _ JMP {k') 

5. If the DPCH spreading factor M < 256 , the above DPDCH channel estimate 
c^uapW is to be interpolated at the DPCH symbol rate to be properly used. The 
least complexity interpolation technique consists of repeating the estimate 
c^-wW 256/Af times, if M £256 and averaging over two successive estimates 
c^uapVc) If M = 512 . In a general preferred embodiment more sophisticated 
interpolation techniques can be adopted (e.g., linear, quadratic, spline interpolation) 

6. Repeat steps 1 to 6 for all paths />= V-.P assigned to the RAKE receiver fingers. 

4.4.2 Joint CPICH-DPCH Low Complexity DPDCH Channel Estimation Algorithm 
In this section we summarize the main steps of the low complexity joint CPICH-DPCH ML 
optimal DPDCH channel estimation algorithm. The following steps are to be implemented to 
estimate the DPDCH channel c^(ft)at time k within the UE. 

1. Compute the ML instantaneous CPICH and DPCH channel estimates c^n) and 
2<m(0 according to equations (8)-(10) and the generalization (23). 

2. Perform interpolation of the above obtained ML Instantaneous DPCH and CPICH 
channel estimates to the lowest symbol rate of said DPCH and CPICH, I.e. 

. If the DPCH spreading factor M< 256 the DPCH ML channel estimates are to be 
interpolated in order to match the CPICH symbol rate 
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• If the DPCH spreading factor M=5l2the CPICH ML channel estimates are to be 
Interpolated in order to match the DPCH symbol rate 

• If the DPCH spreading factor M 256 no Interpolation is needed 

3. The so interpolated ML Instantaneous channel estimates are collected in the following 
vectors 

Kpich 82 \Pcplch (" - Kcplch + ... (riff 

8^*^(11-^+1) ... c dpat (riff 
with / > 0 . In our implementation we shall set / =land K cpkh = K dpch . To further 
reduce the complexity we limit the estimates c^in-K^ +!),» to the 
estimates over the DPCCH c^in-K^ +l),...,c^(«)in order to avoid the 
Implementation of a DD mechanism necessary to provide pilot symbol estimates over 
the DPDCH 

4. Build CPICH auto-correlatlon and DPCH CPICH cross correlation at lag / > 0 of the 
form (17) and (18), i.e., 



and 



5. Build a ML estimate of the beamforming complex factor p as 



h rr\ - $&1D. - 



where we shall set / = 1 according to equation (23). 

6- Build an accurate (e.g. MAP) CPICH channel estimate (k) of the true CPICH 
channel c^Qc) by suitable low-pass filtering the ML CPICH instantaneous channel 
estimates -[^ {n-K^ -Z +1) ... c^{n)f . To reduce the complexity of 
the estimation filter we may consider only a part of the vector of the ML CPICH 
instantaneous channel estimates for low pass filtering. 

7. Build a DPDCH channel estimate c^OQcf the true DPDCH channel as 
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8. If the DPCH spreading factor M <2S6, the above DPDCH channel estimate 

is to be interpolated at the DPCH symbol rate to be properly 
used. As previously mentioned, the least complexity interpolation technique consists 
of repeating the estimate 256/Mdmes, if a/_ 25 6 and 

averaging over two successive estimates ^(*)-A.^(*)lf M=5\2. In a 
general preferred embodiment more sophisticated Interpolation techniques can be 
adopted (e.g., linear, quadratic, spline interpolation) 

9. Repeat steps 1 to 8 for all paths p - I,..., p assigned to the RAKE receiver fingers 

In a non-limitative preferred embodiment, the CPICH multlpath coefficient's estimate at step 
6 is done by a linear prediction as described, e.g„ In the document hereby Incorporated by 
referenced Bsltersee, G. Fock, P. Schultz-RitHch, and H. Meyr, "Performance analysis of 
phasor estimation algorithms for FDD-UMTS RAKE receiver," IEEE 6^ Symp. on Spread 
Spectrum Techonologles and Applications, MJIT, NJ, September 2000." 



4.4.3 Relation to DPCCH Based DPDCH Channel Estimation 

Some authors (see e.g. the hereby Incorporated references "TSG RAN WG 4 meeting #17 
R4-010594 Ericsson, Dedicated Pilots, May 2001," "TSG RAN WG 4 meeting #17, R4-010528 
Nokia, Proposal for user-spedffc beamforming for UTRA FDD, May 2001," and references 
therein, available at http://wvvw.3gpp.org/ftp/tsg_ran/WG4_Radio/TSGR4_17) have 
addressed the problem of DPDCH channel estimation in the presence of transmit 
beamforming. Mainly they conclude that since only the DPCH is affected by the transmit 
beamforming only the Information available from the DPCH channel is to be used for DPDCH 
channel estimation under those circumstances. In particular when DD DPDCH ML channel 
estimates cannot be used due to complexity limitations, they recommend performing DPDCH 
channel estimation only exploiting the a-priori knowledge of the DPCCH pilot symbols. We 
shall notice that when it is viable both the algorithms described here by making use of 
information available from both the CPICH and the DPCCH have two major advantages with 
respect to a DPDCH channel estimation only based on the DPCCH pilots. First, the scarceness 
of pilot symbols provided by the DPCCH (from 1 to a maximum of 16 depending on the slot 
format see the hereby Incorporated reference "3GPP Technical Specifications, Rel. 99, 
http://www.3qpp.org/specs/sp prs html" ), limits the accuracy of DPDCH channel estimate 
because of lack of sufficient noise and Interference suppression. Secondly, as mentioned 
above, within a DPCH time-slot the DPCCH pilot symbols are not transmitted contlnuously as 
the CPICH pilot symbols, but they are time-multiplexed with the unknown data symbols 
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comprising the DPDCH. In the presence of high Doppler spread, I.e. at high UE velocity, such 
structure of the DPCH time-slot might pose the problem of estimating the channel path 
coefficient c dpdch (k) during the absence of known pilot symbols as during the DPDCH period, 

unless the CPICH pilot symbols which are continuously transmitted are not exploited as well. 
5 Indeed If the Doppler spread is very low compared with the slot rate, the channel (*) will 

vary so slowly within a slot period that It could be considered approximately constant over 
that period. Thus the channel estimate provided over the DPCCH can be considered valid for 
the whole DPDCH period. In other words there is no practical need of up-dating the DPDCH 
channel estimate at a higher rate than the slot rate. Conversely if the Doppler spread is of the 

10 same order of magnitude as the slot rate or even larger the assumption that the DPDCH 

channel is approximately constant over a slot period is no longer valid. For Instance referring 
to the UMTS standard, a slot comprises 2560 chip periods yielding a slot rate of 1.5kHz while, 
under extreme cases, the system is supposed to cope with Doppler spread up to 1kHz. 
Therefore in such cases the channel estimated over the DPCCH cannot be considered a valid 

15 approximation of the true channel c^(k) for the entire successive DPDCH period. Under 

these circumstances one should resort to interpolation and prediction techniques or yet 
implement a DD mechanism similar to the one previously described that uses hard decision 
taken on the DPDCH data symbols as known pilots to update the channel estimates In order 
to update at the appropriate rate the channel estimate over the DPDCH period between two 

20 consecutive DPCCH bursts of known pilot symbols. A DPDCH channel estimator, Including 

such a tracking mechanism for the channel estimates is likely too computationally demanding 
to be feasible. Furthermore under certain practical circumstances, such as when the DPCCH 
logical channel consists of only a single pilot symbol per slot, that approach would definitely 
yield poor performances. The joint CPICH-DPCH MAP DPDCH channel estimation algorithm 

25 represents a solution to this problem since it implicitly performs prediction of the true DPDCH 
channel based on the whole available Information adapting the prediction filter to the actual 
Doppler spread and number of DPCCH known pilot symbols. Similarly also the second low 
complexity algorithm proposed here Inherently exploits the presence of the known pilot 
symbols over both the DPCCH and CPICH logical channels and allows tracking of the DPDCH 

30 channel variations even at high Doppler spread. The latter algorithm results more attractive 
from a practical implementation standpoint though, due to the low complexity involved and 
the nearly optimal performance. Indeed the second DPDCH channel estimation algorithm, 

even in its least complexity instance, builds an estimate p of the beamforming complex 
factor p, jointly accounting for the DPCCH and CPICH available information (i.e. all available 
35 known pilot symbols), and provides up-to-date estimates of the DPDCH channel by 

multiplying the estimated factor p by the channel estimate c^(«) computed relying on the 
CPICH pilot symbols. The continuous provision of known pilot symbols over the CPICH 
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channel can be easily exploited by the UE to track the variations due to the Doppler spread of 

the channel path coefficient^) by continuously updating the CPICH channel 

estimate c cpich (k) . The tracking of the DPDCH channel Is inherently obtained as the product 

CfrtchW) *= PttfitptckW) . Finally we shall remark that within our second algorithm the DPCCH 
5 pilots are not directly used to compute a DPDCH channel estimate but only to compute an 
estimate of the factor p. If necessary and/or if the UE can afford the Increased complexity 
hard decision taken on the DPDCH can be used as pilot symbols to Improve the estimate of fi 
as represented by the equations derived above. 

10 4.4.4 Algorithms Working Limitations 

In order to set-up a communication in any mobile communication system, a mobile terminal 
UE at power-up needs to listen to the network and acquire a certain number of parameters 
related to the operator, the BS with which the UE will start the communication, the type of 
service required etc. Prior the useful data transmission the UE has primarily to perform a 
15 certain synchronization procedures to acquire a BS. Within the general framework of DS- 

CDMA systems suitable for implementing a RAKE as receiver at the mobile terminal UE, there 
exists a preliminary synchronization phase that provides along with other parameters and 
time references estimates, also estimates of the time-delays r p , and the average energies 

associated with the paths comprising the channel Impulse response from the BS of interest 

20 In the UMTS standard the mobile terminal UE generally estimates both paths' delays and 
average energies by exploiting the CPICH. Based on those estimates the RAKE receiver 
fingers can be properly set-up, as it Is well known by the person skilled in the art. According 
to the methods of the invention described above, we have Implicitly assumed that when the 
receiver can detect a path and estimate the associated delay on the CPICH during the 

25 synchronization mechanism, a corresponding path at the same delay exists on the DPCH even 
though with a different complex amplitude. As we have seen, this assumption allowed the 
derivation of a DPDCH channel estimation algorithm with low complexity, good performances, 
and high flexibility, able to cope with the problem of DPDCH channel estimation in the 
presence of transmit beamforming. 

30 In the light of the above It will be noted that under certain circumstances the energy 

associated with a path might be high over the CPICH channel and fairly low over the DPCH 
channel. In particular a path can be strong enough to be detected over the CPICH channel by 
the synchronization mechanism, while being too weak to be detected over the DPCH channel. 
Conversely, it may happen that a path that carries significant useful signal energy on the 

35 DPCH channel Is actually undetectable on the CPICH channel, and therefore not detectable by 
the synchronization mechanism operating based only on the CPICH channel. The first case 
would lead to an estimate of p with a very small magnitude. Thus, the algorithms described 




above would still work properly and the RAKE receiver would simply drop the finger as the 
corresponding DPCH channel path's energy becomes too small. The second case is much 
more critical because there would be no way to exploit the Information provided by 
synchronization mechanism operating based only on the CPICH channel. Under these 
5 circumstances any of the two DPDCH channel estimation algorithms would properly operate 
and one should first search for paths associated with the DPCH channel by implementing a 
synchronization mechanism based on the DPCCH. Such a further synchronization mechanism 
would Involve a significant complexity increase along with several real-time management 
complications, and would likely lead to poor performances due to the limited number of 
10 DPCCH pilot symbols per slot as we have previously seen. 

5 DPDCH Channel Estimator Architecture 

In this section we provide architecture schematics for the implementation of the previously 
described DPDCH channel estimation algorithms. 

5.1 Optimal Joint DPCH-CPICH MAP DPDCH Channel Estimation Algorithm 
15 Architecture 

Fig. 3 represents a schematic architecture for an Implementation of a RAKE receiver 
comprising the proposed DPDCH channel estimator using the algorithm summarized by steps 
1—5 described in section 4.4.1. In particular Rg. 3 describes a portion of a single finger of 
the RAKE receiver, comprising: 
20 • A delay recovering block denoted DELAY_RECOVER compensating the delay associated 
with the path assigned to the finger of interest 

• A descrambllng block denoted DESCRAMB 

• A CPICH despreading block, associated with the CPICH channel denoted 
CPICH_DESPREAD 

25 • A DPCH despreading block associated with the DPCH channel denoted DPCH_DESPREAD 

• A CPICH symbol denotation block, denoted CPICH_DEROT, removing the modulation by 
the «-th CPICH pilot symbol 

• A demultiplexer block denoted DEMUX to recover according to the current DPCH slot 
format, the pilot symbols comprising the DPCCH and separate the data symbols 

30 comprising the DPDCH to be decoded 

• A DPCCH symbol denotation block, denoted DPCCH_DEROT, removing the modulation by 
the Mh DPCCH pilot symbol 

• Two interpolation blocks, one for the CPICH denoted INTERP1 and another for the DPCH 
denoted INTERP2. The first Interpolation block INTERP1 Is enabled when the DPCH 

35 spreading factor M = 5 12 . The second Interpolation block INTERP2 is enabled only when 
the DPCH spreading factor M < 256 . No interpolation block is enabled when M = 256 (see 
step 2 section 4.4.1). 
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• A joint CPICH-DPCH MAP DPDCH channel estimation low pass /7/ier denoted 
DPDCH_MAP_CH_EST, involving a delay of L CPICH symbol periods, and providing the 
estimates c^-** (*) according to step 4 of section 4.4.1. 

. A delay line of L(256/iif) DPCH symbol periods denoted DELAY necessary to temporally 
align the MAP DPDCH channel estimate c^.^, (ft) with the DPDCH data symbols 

. An interpolation block denoted INTERP3 to adapt the estimate at the CPICH symbol rate 
to the DPCH symbol rate, according to step 5 of section 4.4.1. 

• A block denoted CMPLX_CNJ, computing the complex conjugate of the MAP DPDCH 
channel estimate after interpolation to the DPCH symbol rate at the output of the block 
INTERP3, necessary to the RAKE combiner, as It is well known by tiie person skilled in the 
art 

. A multiplier block denoted DPDCH_CHEST_DATA_MULT that computes the product of the 
time-aligned DPDCH symbols and the corresponding complex conjugate DPDCH channel 
estimate. The product result is then sent to the RAKE receiver combiner, denoted 
RAKEjCOMBINER. 

The scheme depicted in Rg. 3 it to be replicated for all RAKE fingers. The output of the 
multiplier block DPDCH_CHEST_DATA_MULT of each finger Is summed up by the RAKE 
combiner RAKEjCOMBINER to form the output of the RAKE receiver as It Is well known by 

the ordinary skilled person. 

In practice, In a preferred embodiment, nearest neighbor Interpolation will be adopted for 
the implementation of all the interpolation blocks INTERP1, INTERP2, and INTERP3, as it is 
shown to be a proper solution for this purpose In terms of complexity and performance trade- 
off. Nearest neighbor interpolation Is the least complexity interpolation method. It actually 
saves complexity at expense of a higher signal distortion than that resulting with more 
sophlsticat6ed Interpolation techniques. However it can be shown that the actual amount of 
signal distortion yielded by that interpolation method is significantly below the noise and 
interference for practical ranges of signal, interference, and noise levels. Therefore the 
additional signal distortion due to nearest neighbor interpolation is actually negligible 
compared with the effects of the noise and interference so as not affecting the receiver 
performances. 

5.2 Joint CPICH-DPCH Low Complexity DPDCH Channel Estimation Algorithm 
Architecture 

Rg. 4 represents a schematic architecture for an Implementation of a RAKE receiver 
comprising the proposed DPDCH channel estimator using the algorithm summarized by steps 
1—8 described in section 4.4.2. In particular Rg. 4 describes a portion of a single finger oi 
the RAKE receiver, comprising: 
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A delay recovering block denoted DELAY.RECQVER compensating the delay associated 
with the path assigned to the finger of Interest 
A descrambllng block denoted DESCRAMB 

A CPICH despreadlng block, associated with the CPICH channel denoted 
CPICH_DESPREAD 

A DPCH despreadlng block associated with the DPCH channel denoted DPCH_DESPREAD 
A CPICH symbol derotation block, denoted CPICH.DEROT, removing the modulation by 
the n-th CPICH pilot symbol 

A demuldplexerUock denoted DEMUX to recover according to the current DPCH slot 
format, the pilot symbols comprising the DPCCH and separate the data symbols 
comprising the DPDCH to be decoded 

A DPCCH symbol derotation block, denoted DPCCHJ5EROT, removing the modulation by 
the Hh DPCCH pilot symbol 

A CPICH channel estimator low pass /7/terdenoted CPICH_CH_EST, Involving a delay of 
L CPICH symbol periods, and providing the estimates c^Qc) according to step 6 of 
section 4.4.2 

The actual estimator of the beamforming factor ft, which comprises two interpolation 
blocks, one for the CPICH denoted INTERP1 and another for the DPCH denoted INTERP2. 
The first interpolation block INTERP1 is enabled when the DPCH spreading factor 
M m 512 . The second Interpolation block INTERP2 Is enabled only when tine DPCH 
spreading factor M < 256 . No Interpolation block Is enabled when M = 256 see step 2 
section 4.4.2. 

. A delay line of L CPICH symbol periods denoted DELAY1 necessary to temporally align 
the CPICH channel estimate c^(*)with the up-to-date estimate ft of the beamforming 
factor ft 

• A delay line of i(256/M) DPCH symbol periods denoted DELAY2 necessary to temporally 
align the CPICH channel estimate (ft) with the DPDCH data symbols 

• A low-pass-fllter denoted LPF1 for the DPCH-CPICH cross-correlation estimate and a LPF 
denoted LPF2 for the CPICH autocorrelation estimate to be enabled when the UE Is aware 
that the beamforming factor ft stays constant The LPF's memories are to be reset as the 
UE detects that the beamforming complex factor has abruptly changed (either via 
signaling from the BS or via a suitable detection mechanism implemented within the UE). 
Both filters LPF1 and LPF2 are not to be used if the UE has no mean to detect abrupt 
changes of the beamforming complex factor. 

• A multiplier denoted BETA_CPICH_CH_EST_MULT computing the product fi^c^ (k) . 
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• An interpolation block denoted INTERP3 to adapt the estimate at the CPICH symbol rate 
? **s(*) = Puic^tk) obtained at step 7 to the DPCH symbol rate, according to step 8 of 
section 4.4.2. 

• A block denoted CMPIX.CNJ, computing the complex conjugate of the DPDCH channel 
estimate after interpolation to the DPCH symbol rate at the output of the block INTERP3, 
necessary to the RAKE combiner, as it is well known by the person skilled in the art 

• A multiplier block denoted DPDCH_CHEST_DATA_MULT that computes the product of the 
time-aligned DPDCH symbols and the corresponding complex conjugate DPDCH channel 
estimate. The product result Is then sent to the RAKE receiver combiner, denoted 
RAKE_COMBINER. 

The scheme depicted in Rg. 4 it to be replicated for all RAKE fingers. The output of the 
multiplier block DPDCH_CHEST_DATA_MULT of each finger is summed up by the RAKE 
combiner RAKE_COMBINER to form the output of the RAKE receiver as It is well known by 
the ordinary skilled person. 

In practice, in a preferred embodiment, also in this case nearest neighbor interpolation 
will be adopted for the implementation of all the interpolation blocks INTERP1, INTERP2, and 
INTERP3, as it is shown to be a proper solution for this purpose in terms of complexity and 
performance trade-off. Indeed the same considerations on nearest neighbor interpolation 
done in section 5.1 on the joint CPICH-DPCH MAP DPDCH channel estimation algorithm hold 
in the case of the joint CPICH-DPCH low complexity DPDCH channel estimation algorithm as 
well. 

6 Concluding Remarks 

It is to be understood that the present invention is not limited to the aforementioned 
embodiments and variations and modifications may be made without departing from the spirit 
and scope of the invention as defined in the appended dalms. In this respect, the following 
concluding remarks are made. 

It Is to be understood that the present Invention is not limited to the aforementioned 
UMTS application. It can be used within any application of DS-CDMA where the 
communication standard envisages the presence of two logical channels, i.e. a common pilot 
channel, comprising symbols or signal waveform(s) known at the mobile terminal UE, and a 
dedicated channel comprising unknown data symbols and, In a preferred but not limitative 
embodiment, also dedicated known pilot symbols or signal waveform(s). 

It has to be understood that the present invention can be applied in the uplink 
communication, from the mobile terminal UE to the base station BS, when the mobile 
terminal UE the employs transmit beamforming. 




It is to be understood that the methods according to the present Invention are not 
limited to the aforementioned Implementation. There are numerous ways of Implementing 
functions of the methods according to the invention by means of items of hardware or 
software, or both, provided that a single item of hardware or software can carries out several 
5 functions. It does not exclude that an assembly of items of hardware or software or both 
carry out a function. For instance, algorithms steps can be combined, thus forming a single 
function without modifying the methods of channel estimate with transmits beamforming in 
accordance with the invention. 

Said hardware or software items can be implemented in several manners, such as by 
10 means of wired electronic circuits or by means of an Integrated circuit that Is suitably 
programmed respectively. The integrated circuit Is to be part of a receiver that can be 
contained in a computer, in a mobile communication handset, in a base-station, or In any 
other communication system device. Said receiver comprises means of the receiver adapted 
to make all necessary operations to support a certain type of communication services (I.e. 
15 operations of de-scrambling, de-spreading, etc.) said means being hardware or software 
items as above stated. In a similar fashion, the first estimator within the receiver comprises 
maximum likelihood estimation means adapted to make the ML Instantaneous estimates of 
step 1 and 3, interpolation means adapted to make rate matching as in step 2, optimal MAP 
DPDCH channel multipath coefficient estimation means adapted to make the DPDCH 
20 estimates of step 4, interpolation means adapted to make rate matching as in step 5, as 
previously described, said means being hardware or software items as above stated. In the 
same way, the second estimator within the receiver comprises maximum likelihood estimation 
means adapted to make the ML instantaneous estimates of step 1 and 3, interpolation means 
adapted to make rate matching as in step 2, correlation estimation means adapted to make 
25 the correlations estimates of step 4, beamforming complex factor estimation means that are 
adapted to make the factor complex estimate of step 5, a first multipath coefficient 
estimation means adapted to make the said CPICH multipath coefficient estimate of step 6, a 
second multipath coefficient estimation means that are adapted to make said DPDCH 
multipath coefficient estimate of step 7, Interpolation means adapted to make rate matching 
30 as In step 8, as described previously, said means being hardware or software items as above 
stated. Of course, as mentioned above, said means can be an assembly of items of items of 
hardware or software or both carry out a function, or a single Item carrying out several 
functions. 

The Integrated circuit comprises a set of Instructions. Thus, said set of instructions 
35 contained, for Instance, In a computer programming memory or in a stand-alone receiver 
memory may cause the computer or the receiver to carry out the different steps of the 
estimate methods. 
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The set of instructions may be loaded into the programming memory by reading a 
data carrier such as, for instance, a disk. A service provider can also make the set of 
instructions available via a communication network such as, for example, the Internet. 

Any reference sign in the following claims should not be construed as limiting the 
claim. It will be obvious that the use of the verb "to comprise" and its conjugations do not 
exclude the presence of any other steps or elements besides those defined in any claim. The 
article "a" or w an" preceding an element or step does not exclude the presence of a plurality 
of such elements or steps. 



CLAIMS 

Method for estimating a propagation channel In the presence of transmit beamforming, 
accounting for the structure of two logical channels (CPICH, DPCH) and based on a 
common structure of corresponding propagation channels, said second logical channel 
(DPCH) comprising two sub-channels (DPDCH, DPCCH), said propagation channels being 
modeled as a linear superposition of a finite number of discrete multipath components 
(p=l,...,P) following an uncorrelated-scattering wide-sense stationary model, a multipath 
component being characterized by a time-varying multipath complex coefficient 
( Cp 0) and fi p c p (t) ) and a delay ( r p )• 

A method for estimating propagation channel in the presence of transmit beamforming as 
claimed in dalm 1, characterized in that said propagation channel correspond to the first 
sub-channel (DPDCH) and that said method provides estimates of each multipath 
component (p=l>~., P) complex coefficient ( P p c p (t) ) according to a max!mum-a- 
posteriori MAP optimization criterion accounting for the whole available information 
associated with said logical (CPICH, DPCH) and corresponding propagation channels, 
through the following processing steps: 

1. Build a second channel (DPCH) and a first channel (CPICH) Instantaneous 
maximum likelihood ML channel multipath complex coefficients estimates 

(»)), 

2. Perform Interpolation of the above obtained ML instantaneous second (DPCH) 
and first (CPICH) channel multipath complex coefficient estimates (£4*61) and 
c^in) ) to the lowest symbol rate of said second (DPCH) and first (CPICH) 
logical channels, 

3. Compute an optimal linear prediction filter (f ) according to a joint second and 
first channels (DPCH-CPICH) maxlmum-a-postertorl (MAP) criterion, 

4. Filter the Interpolated ML Instantaneous second (DPCH) and first (CPICH) channel 
multipath complex coefficient estimates obtained at step 2 with said optimal 
linear prediction filter in order to obtain a MAP first sub-channel (DPDCH) 
multipath coefficient estimate {.c JpJch . UAP {k)) l and 

5. Interpolate said MAP first sub-channel (DPDCH) multipath coefficient estimate 
(e*rf*-*«p(*>)to the second logical channel (DPCH) symbol rate when said 
symbol rate Is lower than the first logical channel (CPICH) symbol rate. 

where steps 1 to 5 are repeated for ail multipath component (p^l,...?) complex 
coefficients (/? p c p (0). 
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3. A second method for estimating a propagation channel In the presence of transmit 
beamformlng characterized in that said propagation channel corresponds to the first sub- 
channel (DPDCH) and that said method provides estimates of each multipath component 
(p=l,...,P) complex coefficient, accounting for the whole available information associated 
with said logical (CPICH, DPCH) and corresponding propagation channels, through the 
following processing steps: 

1. Build a second channel (DPOI) and a first channel (CPICH) instantaneous 
maximum likelihood ML channel multipath coefficients estimates (c^^and 

2. Perform interpolation of said ML instantaneous first (DPCH) and second (CPICH) 
channel multipath coefficient estimates (c^OOand c^Qi) ) to the lowest 

symbol rate of said second (DPCH) and first (CPICH) logical channels, 

3. Build an optimal maximum a posteriori estimate (c^-awpW) of ^ e flrst 
(CPICH) channel multipath coefficient {c^ teh (fc) ), 

4. Build an estimate of the cross-correlation ( $ dc (/) ) between the first (CPICH) and 
second (DPCH) channel multipath coefficient Instantaneous maximum likelihood 
estimates obtained at step 2 {c dpch and ) and an estimate of the 

autocorrelation ( (I) ) between the (CPICH) channel multipath coefficient 
instantaneous maximum likelihood estimates (c^ /dk ) of step 1 and 2 at non-zero 
correlation lag ( / * 0 ) for noise suppression, 

5. Build an estimate ( p ) of a beamformlng complex factor ( p ) from said cross- 
correlation and autocorrelation estimates, 

6. Build a first sub-channel (DPDCH) multipath coefficient estimate ( c dpdch (k) ) as 

the product of the estimates obtained at steps 3 ( c^ MAP (k) ) and 5 ( p ), and 

7. Interpolate said first sub-channel (DPDCH) multipath coefficient estimate 
(c dpdch (k) ) to the second logical channel (DPCH) symbol rate when said symbol 

rate Is lower than the first logical channel (CPICH) symbol rate 

where steps 1 to 7 are repeated for all multipath component (p=l,...,P) complex 
coefficients (fi p c p (Q). 

4. A method as claimed In claims 2 and 3, characterized In that the first logical channel 
(CPICH) maximum likelihood channel multipath coefficient estimates (c^fa) ), are 
computed based on the a-prlori knowledge of some pilot symbols forming said first logical 
channel (CPICH). 




5. A method as claimed in claims 2 and 3, characterized in that the second logical channel 
(DPCH) maximum likelihood channel multipath coefficient estimates {c dpch {n) ), related to 

the second sub-channel (DPCCH) are computed based on the a-priori knowledge of the 
pilot symbols forming said second sub-channel (DPCCH). 

5 6. A method as claimed in claims 2 and 3, characterized In that the second logical channel 
(DPCH) maximum likelihood channel multipath coefficient estimates (c^fr)), related to 
the first sub-channel (DPDCH) are computed by a decision-direct mechanism. 

7. A method as claimed in claims 2 and 3, characterized In that the interpolation of step 2 is 
performed by nearest neighbor interpolation. 

10 8, A method as claimed in claim 2, characterized In that the optimal linear prediction filter Is 
built according to the maximum-a-posteriori optimization criterion, based on the 
interpolated maximum likelihood channel multipath coefficients estimates (c^O*) and 

^(n)) related to said first (CPICH) and second (DPCH) logical channels in order to 

provide an optimal joint second and first channels (DPCH-CPICH) maximum-a-posteriorl 
15 first sub-channel (DPDCH) multipath coefficient estimate ( c^^^p (k) ). 

9. A method as claimed in daim 3, characterized In that a maximum likelihood estimate of 
the second (DPCH) corresponding propagation channel and first (CPICH) corresponding 
propagation channel cross-correlation (B{c dpeh GOc^ (*-/)}) and a maximum likelihood 
estimate of the first (CPICH) corresponding propagation channel autocorrelation 

20 ( E^,^ (7i)c^ (it - 1)} ) are computed based on the sample moments ( 4.(0 and 4, (/) ) 
of the first (CPICH) and second (DPCH) channel maximum likelihood estimates 
(£«ipiOO and c dpch {n)) of step 1 and 2. 

10. A method as daimed in daim 3, for the computation of the estimate of said complex 
beamforming factor (fi) characterized in that the second logical channel (DPCH) and the 

25 first logical channel (CPICH) corresponding propagation channel cross-correlation and the 
first logical channel (CPICH) corresponding propagation channel autocorrelation 
maximum likelihood estimates ( f+ (/) and ^ (/) ) at different correlation lags (£=l, 2,. . .JL) 

are linearly combined ( a,fa{l) and ^ b£ m (0 )• 



36 



11. A method as claimed in claim 3, characterized In that the second logical channel (DPCH) 
and first logical channel (CPICH) cross-correlation and the first logical channel (CPICH) 
autocorrelation successive estimates (&,(/) and^(/» are taken at a fixed lag (0 and 
are low-pass filtered for the computation of the estimate of said complex factor ( p ). 

5 12. A method as claimed in claim 3, characterized in that the estimate of said complex factor 
( p ) is built as a linear combination of the beamforming complex factor estimates 
computed as the ratio of the second logical channel (DPCH) and the first logical channel 
(CPICH) corresponding propagation channels cross-correlation and the first logical 
channel (CPICH) corresponding propagation channel autocorrelation estimates at a 

LO certain lag (0 (A*(0 = £>(D/<L(0), (£=£r,A*(0) at,a 9 '=1.2,...,*. 

13. A method as claimed in any claims 10, 11 or 12, characterized in that the estimate of said 
complex factor ( p ) is limited to the lag equal to 1. 

14. a receiver utilizing said methods as claimed in any claims 1, 2 or 3. 

An Estimator for estimating a propagation channel in the presence of transmit 
beamforming, accounting for the structure of two logical channels referred as to a 
common channel and a dedicated physical channel (CPICH, DPCH) and based on a 
common structure of corresponding propagation channels, said dedicated physical 
channel (DPCH) channel comprising two sub-channels (DPDCH, DPCCH), said 
propagation channels being modeled as a linear superposition of a finite number 
(p=i„..,p) of discrete multlpath components following an uncorrelated-scattering wide- 
sense stationary model, a multlpath component being characterized by a time-varying 
multlpath complex coefficient (c,(0 and p p c p (t) ) and a delay (t„). 



16. An estimator as claimed in claim 15 for estimating a propagation channel In the presence 
of transmit beamforming, characterized in that said propagation channel corresponds to 
the first sub-channel (DPDCH) estimation and that said estimator comprises: 
. Means to build a second (DPCH) and a first (CPICH) logical channel corresponding 

propagation channels instantaneous maximum likelihood ML channel multlpath 

coefficients estimates ( (m) and («) ), 
. Means to perform Interpolation of the above obtained ML instantaneous second 

(DPCH) and first (CPICH) logical channel corresponding propagation channels 
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multipath coefficient estimates (c^(n) and a^OO) to the lowest symbol rate of 
said second (DPCH) and first (CPICH) logical channels, 

. Means to build an optimal linear prediction filter according to a Joint second and first 
(DPCH-CPICH) channels maxlmum-a-posterlori criterion, 

. Means to build a first sub-channel (DPDCH) multipath coefficient estimate 
{c ifdMMP Qc)) by filtering with said optimal linear prediction filter with said 
interpolated ML Instantaneous second (DPCH) and first (CPICH) logical channel 
corresponding proagatlon channels multipath coefficient estimates ( and 

),obtalned at step 2, and 
. Means to Interpolate said first sub-channel (DPDCH) multipath coefficient estimate 
(w P W)to the second logical channel (DPCH) symbol rate when said symbol 
rate is lower than the first logical channel (CPICH) symbol rate. 

17. An estimator as claimed In claim 15 for estimating a propagation channel In the presence 
of transmit beamforming, characterized In that said propagation channel corresponds to 
the first-sub-channel (DPDCH) and that said estimator comprises: 

• Means to build a second (DPCH) and a first (CPICH) logical channel corresponding 
propagation channels instantaneous maximum likelihood ML channel multipath 
coefficients estimates ( («) and (n) ), 

• Means to perform interpolation of the above obtained ML instantaneous second 
(DPCH) and first (CPICH) logical channel corresponding propagation channels 
multipath coefficient estimates (c^OO and ^(n)) to the lowest symbol rate of 
said second (DPCH) and first (CPICH) logical channels, 

• Means to build an optimal maximum a posteriori estimate ( c qM ^. IMP (*) ) of the first 

logical channel (CPICH) multipath coefficient (c^ (A)), 

• Means to build an estimate (4(0) of the cross-correlation (E^Wc^ («-/)} ) 
between the first (CPICH) and second (DPCH) logical channel corresponding 
propagation channels multipath coefficients instantaneous maximum likelihood 
estimates (c^OO and c^(n)), and an estimate (4(0) of the autocorrelation 
(E{c <pWi (n)c; Wi (n -/)} ) between the first logical channel (CPICH) corresponding 
propagation channel multipath coefficient instantaneous maximum likelihood 
estimates ( (/») ), of step 1 and 2 of daim 3, at non-zero correlation lag ( / * 0 ) for 

noise suppression, 

• Means to estimate a beamformlng complex factor ( p ) from said cross-correlation 
and the auto-correlation estimates (4(0 and 4(0)/ 
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• Means to build a first sub-channel (DPDCH) multipath coefficient estimate ( c dpdch (k) ) 
as the product of the optimal maximum a posteriori estimate (c tpkl ^ ]li4P (k)) of the 
first channel (CPICH) multipath coefficient and the cross-corelatlon and the auto- 
correlation estimates (fat) andfal)) 
5 • Means to interpolate said first sub-channel (DPDCH) multipath coefficient estimate 
(Ccptcb-MApOO) to the second logical channel (DPCH) symbol rate when said symbol 
rate is lower than the first logical channel (CPICH) symbol rate 

18. A receiver comprising an estimator as claimed in claim 15. 

19. A communication system using a method for estimating a propagation channel In the 
10 presence of transmit beamforming as claimed in claim 1, when information data are 

transmitted through a beamforming system. 
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Methods for channel estimation in the presence of transmit beamformlng 
ABSTRACT 

The present Invention relates to a method for estimating a propagation channel in the 
presence of transmit beamformlng within a receiver, accounting for the structure of two 
5 logical channels referred as to a common channel and a dedicated physical channel (CPICH, 
DPCH) and based on a common structure of corresponding propagation channels, said DPCH 
channel comprising two sub-channels (DPDCH, DPCCH), said propagation channels being 
modeled as a linear superposition of a finite number of discrete multipath components 
following an uncorrelated-scattering wide-sense stationary model, a multipath component 
10 being characterized by a time-varying multipath complex coefficient and a delay. 

Use: RAKE Receiver 

Reference: Fig. 3 
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